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Motion of classical point-like impurities in the homogeneous Einstein condensate of bosons is 
studied in the framework of second quantization method. A toy model is proposed and its general 
solution within the Bogoliubov approximation is obtained. The effective Minkowski space-time 
structure arises naturally in this non-relativistic quantum many-body system in the low energy 
regime. This is shown to be true in this model. Several examples are discussed in order to illustrate 
our model. The homogeneous condensate produces an effective Yukawa type attractive force between 
impurities sitting in condensate. Landau's criterion is naturally derived in a case of linear motion of 
impurity. The analytic expressions for spectra of Bogoliubov excitations produced by the accelerated 
motions of impurities are obtained. A quick look at the analytic expression reveals that the spectrum 
of gapless excitations emitted by the linearly accelerated impurity is not thermal. If the homogeneous 
condensate is the physically correct model for Minkowski space-time then it follows that the apparent 
thermal response of the simple linearly accelerated detector models may be the result of improper 
regularization. 

PACS numbers: 03.75.Kk, 03.70.+k, 67.40.Yv 



I. INTRODUCTION 

In 1924 Einstein predicted the very peculiar behavior of an ideal gas of non-relativistic particles obeying a new 
statistics [1]. This new statistics (Bose-Einstein statistics) was proposed by Bose to explain the Planck distribution for 
the black body radiation. Einstein found that most of the population of the ideal gas would occupy the ground state 
below the critical temperature. This phenomenon is somewhat counterintuitive, and has no place in classical statistics, 
but it revealed the quantum nature of many particle systems. The phenomenon is now known as Bose-Einstein 
condensation, referring to the Einstein condensation (entartung) in systems having the Bose-Einstein statistics. In 
1995 Bose-Einstein condensate (BEC) was finally created in laboratories in several alkali vapors [2-4]. The essential 
difference besides developments in cooling techniques is an additional magnetic trapping to confine the cold atoms 
within a small region. This celebrated experimental success also revived theoretical interests in BECs with additional 
trapping potentials [5]. 

A decade after the experimental realization of BEC, this field is still very active, both experimentally and theo- 
retically. Recently much interest has been growing in cosmology and gravitational physics. In these fields it is very 
difficult to do experiments in laboratories, and experimental tests are usually based on astrophysical observations. 
There are several proposals to resolve this discrepancy, which discuss "analogies" between real gravity and BECs [6] . 
In their proposals the homogeneous BEC serves as a Minkowski space-time and the deviations from the flat metric 
tensor arise as a so called "acoustic metric" in inhomogeneous BECs. In this way BECs may mimic "analog" gravities. 
Hence "effective" gravities in laboratories are attainable to test the theories. 

On the other hand, there is another direction of research where gravitation is indeed considered as an emergent 
phenomenon arising from collective excitations in quantum many-body systems [7-15]. In these proposals a physical 
vacuum is considered as a real condensed state which contains a huge number of heavy bosons. The Minkowski vacuum 
is assumed to be the homogeneous BEC. The low energy vacuum structure of such an medium is then described by 
the relativistic quantum field theory. In this sense the relativistic quantum field theory is an effective theory which 
works only in the low energy scale. It then follows that the relativistic quantum field theory breaks down at certain 
energy scale. Many difficulties of the relativistic quantum field theory seem to occur as a result of an extrapolation 
and the usual regularization techniques corresponding to unphysical degrees of freedom. 

The idea of the space-time medium for the propagation of light in quantum mechanics was introduced by Dirac 
[16, 17] . The old idea of an aether was abandoned as a result of success of special relativity and by the Michelson-Morley 
experiment before quantum mechanics. Recent precise experiment shows that the upper limit for the anisotropy of 
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the velocity of light is 1CP 15 [18]. Therefore it is an experimental fact that the velocity of light is isotropic. Dirac 
revisited the possibility of keeping the idea of an aether in quantum mechanics. He showed that the isotropy of the 
velocity of light does not conflict with the idea of the space-time medium which is responsible for the propagation 
of light. From this point of view, a quantum many-body system could be a candidate since the velocity of sound is 
isotropic in several BECs and in 4 Hc. 

It is of fundamental interest to study the motion of impurities in BECs. By impurities we mean any bosons or 
fermions which can be distinguishable from condensed bosons, and the total number of impurities N imp is much less 
than that of condensed bosons N (N 3> N imp ). The general solution to this problem has not been found so far, and 
we are unlikely to obtain it within our current understanding in many-body problems. Several partial results are 
obtained recently from a condensed matter physics point of view [22-28]. This problem is not only interesting in a 
pure academic sense. It is also very important for practical reasons to take into account the effects of impurities in 
experiments. 

The purpose of this paper is to provide a small step toward solution of this problem. We propose a toy model 
describing the motion of classical impurities in the homogeneous BEC in the framework of the second quantization 
method. To this end we employ the Bogoliubov model for weakly interacting massive bosons in a dilute regime 
[29-31]. We then add an interaction Hamiltonian which takes into account the local contact interaction between 
bosons and impurities. Specifically, we focus on the motion of classical point-like impurities moving along given 
trajectories. We do not discuss back reaction on the impurities, and their dynamics mainly. Our main interest is 
to find out how the homogeneous condensate responds to the motion of impurities, and how it tries to maintain its 
condensed state in the presence of impurities. We also want to examine the hypothesis postulated in previous studies, 
stating that the homogeneous BEC is the Minkowski vacuum. Especially important is finding out how the effective 
Minkowski space-time structure emerges from the non-relativistic many-body system. Is the motion of impurities in 
the homogeneous BEC exactly corresponding to the motion of test particles in the Minkowski space-time? Suppose 
that the homogeneous condensate happens to be unstable against the motion of impurities and can be easily destroyed 
by them. If it happens, then their hypothesis should be rejected. This is because the Minkowski vacuum has to be 
very stable, as is well known. We will show that the homogeneous condensate is very robust in a large N limit and 
can give rise to macroscopic effects which could be observed. 

After we demonstrate that the homogeneous condensate has similar properties to the Minkowski vacuum, we will 
study several examples to illustrate our model. Particularly, the accelerated motion of impurities in BEC is of interest. 
This situation is analogous to the case of accelerated moving detectors in the Minkowski space-time. According to 
the standard description of quantum field theory in the curved space-time, a uniformly accelerated detector in the 
Minkowski space-time would detect a thermal spectrum. In other words, observers "feel hot" when they are accelerated 
in the Minkowski space-time. The temperature T a which such an observer would measure is given by the formula, 

T «=2^- « 

Here a is the acceleration constant of the observer at an instantaneous rest frame, and ks is the Boltzmann constant. 
This thermal effect is generally referred to as the "Unruh effect" [32]. Thus, it has long been believed that an empty 
Minkowski vacuum would serve as a thermal bath for accelerated observers. In spite of the absence of experimental 
confirmation, this mathematical consequence has been accepted for over three decades. There is sometimes an 
impression that this thermal character is equivalent to quantum field theory itself [33]. Recently, Belinskii and his 
collaborators revisited the Unruh effect. They claim that the thermal character of the Minkowski space is an artifact 
of an invalid mathematical procedure repeatedly used in the discussion of the Unruh effect [34-36]. The problem 
remains open. 

Therefore it is very interesting to examine the accelerated motion of impurities in the homogeneous condensate. In 
the model presented in this paper, the thermal spectrum is not found for cases such as a uniform accelerated motion, 
nor a uniform circular motion. This paper is organized as follows. In Sec. II we propose a simple model Hamiltonian 
describing a point-like classical impurity in BEC, and we obtain the general solution within the Bogoliubov approx- 
imation. In order to keep this paper self-contained, we give the review on the number conserving Bogoliubov model 
briefly in this section. We then derive formulae for the occupation number and the energy dissipation in Sec. III. The 
general solution obtained in Sec. II provides time dependent behavior of various physical quantities such as the den- 
sity, the current density, the stress energy momentum complex. Deviations of these quantities from the homogeneous 
condensate are also discussed. Low energy behavior of these quantities suggest that the Minkowski space-time be an 
emergent space-time arising from non-relativistic quantum many-body system, i.e. the homogeneous BEC. In Sec. IV 
we study several examples to illustrate our model and its general solution. We first provide the general argument for 
the derivation of Landau's criterion based on the method of stationary phase. Next we obtain an attractive Yukawa 
type force between two stationary impurities sitting in the BEC. Then, three particular examples for the motion of 
impurity are studied; a linear motion with a constant velocity, a uniformly accelerated motion, and a circular motion 
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with a constant angular velocity. We obtained analytic expressions for these cases with thorough analyses of them. 
Sec. V provides brief discussion for the dynamics of a quantum impurity in our model. We obtained a closed equation 
for the motion of single, point-like quantum impurity. In Sec. VI we also give an alternative model for the motion 
of impurities based on the Galilean invariance. The current-current type interaction model is proposed and solved 
within the same approximation done in Sec. II. We summarize our results and give a discussion on our results in 
Sec. VII. Appendix contains a supplement for the matin text. 



II. MODEL HAMILTONIAN AND ITS DIAGONALIZATION 

A. The Model 

The Bogoliubov model Hamiltonian for weakly interacting bosons of mass M is 

H (t) = Jd 3 x^(x,t)(-^ f )^(x,t) + ^ jd 3 xjd 3 x , ^(x,t)^(x , ,t)V(x-x')^(x,t)^(x',t). (2) 

The field operators ffi(x,t) and 4>(x,t) in the Heisenberg picture satisfy the equal time canonical commutation 
relations : 

$(x,t),ft(x',t)]=5(x-x'), 
$(x,t),$(x',t)} = [ft(x,t),ft(x',t)] = 0. 

The general two-body interaction term V(x — x') can be approximated by a contact interaction gS(x — x') in the low 
energy region. This approximation will be used throughout the paper in order to get analytic results, but extension 
to the general two-body interaction is quite straightforward [37]. The coupling constant g for the approximated 
interaction between bosons is expressed in terms of the s-wave scattering length a s ; g = Aita s h 2 /M . We assume the 
repulsive interaction g > 0. The condition of diluteness na 3 <gc 1 is also assumed. Here n is the number density of 
massive bosons. We consider the following interaction Hamiltonian between bosons and the distribution of classical 
impurities p c (x, t) taking into account the local interactions between them. 

H\{t) = X J d 3 x p c {x,t)ft{ x,t)ip(x,t). (4) 

For a point-like classical impurity moving along a given trajectory p c (x,t) is written using the Dirac delta 
distribution as p c (x,t) = 5{x — C{t)). The coupling constant A is expressed in terms of the s-wave scattering length 
b s ; A = 2irbl h 2 /m. Here m is the reduced mass of the boson and the impurity. Extension to an arbitrary number of 
point-like impurities case is obtained by a replacement Xp c (x,t) — > J2 i XiS(x — Ci{t))- 

The momentum operator of the bosons 

P = J d 3 x ft(x,t)(-ihV)ip(x,t), (5) 

commutes with the Hamiltonian (2). However, the momentum operator does not commute with the interaction Hamil- 
tonian (4). This is contrary to the fundamental requirement of translational invariance. This should be understood 
as follows. Let us consider a quantum mechanical impurity described by its position and momentum operator q(t) 
and p(t) respectively. The interaction Hamiltonian in this case reads 

Hi(t)=X J d 3 x5(x-q(t))ft(x,t)i>(x,t). (6) 

Then the total momentum operator P + p commutes with this interaction Hamiltonian. The non-conserving mo- 
mentum is now understood as an approximation of replacing q-number with c-number impurity trajectory £(t). This 
approximation is valid as long as only slowly moving classical impurities are concerned. The quantum impurity will 
be discussed later in Sec. VI, where the closed equation is obtained for the quantum impurity. 

We expand the field operators in terms of the plane wave basis with periodic boundary conditions in a finite size 
box V = L 3 as 

Ho(t) =J2^a{(t)a k (t) + 4+,(*)4'-,(*)6fc(*)M*) (7) 

k k,k' ,q 
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and 

H,(t\ = 

V 



#i(*) = ^I>*'(*)4+*'(*) a *(t), ( 8 ) 

k,k' 



where = h 2 k 2 /2M is the free kinetic energy of bosons. Here the summation is taken over integers n, for k = 
2ir(n x , n y , n z )/L. In the above expression, we introduced the Fourier transform of the impurity distribution : 

Pk(t) = J d 3 x p c (x,t)e- ikx . (9) 

By definition p~k=o(t) = J d 3 xp c (x,t) = 1. 

The creation and the annihilation operators a k (t), dk(t) satisfy the equal time canonical commutation relations : 

[d k (t),a{,(t)}=5 k ,k>, ^ 
[fi*(t), &*'(*)] - [a{(t),al(t)] = 0, 

where 5k,k' — 5 nx ,n' &n y ,n' $n z ,ri with <5„ >n / the Kronecker delta. In the following discussion, time dependence of these 
operators for simplicity will be suppressed. 

Later we will take the thermodynamic limit in which the volume V and the number of bosons N become infinite, 
keeping the density n = N/V fixed. In this limit the summation is replaced with the integral 

^E = ( 2 -r 3 / d3fc - (") 

The creation and annihilation operators are relabeled by the continuous wave numbers as 

lim Wat = (2 7 r) 3 / 2 a t (fe), lim Wa k = (2n) 3 ^d(k). (12) 



Then, a^(fe) and a(k) satisfy the canonical commutation relations: 

[a(k), tf(k')] =6{k-k'), 
[a(k), a(k')} = [tf(k), a*(k')] = 0. 



(13) 



In the following we are interested in the effects of impurities to the homogeneous condensate in the large N limit. Wc 
will investigate the first order correction on the homogeneous condensate. 

We follow Bogoliubov's treatment [29-31] to simplify the full Hamiltonian Hq + H\. First we separate the second 
term in (7) into two parts, one is the forward scattering term with zero momentum transferred J7o and the other with 
non zero momentum transferred U' . Uo and U' are 

^o = ^E44«fe«fe' = £(# 2 -#)> ( 14 ) 

k,k' 

where N = ^ fe a k dk is the total number operator, and 

V' = W E 4+<,4'-<, a fe«fe'- ( 15 ) 

fc,fc',q#0 

The second term U' is split into three parts depending on the number of the zero mode creation and annihilation 
operators as U' = + U<V + jjM : 

U {2) = ^al&o + ^7( a o«o 4 a -fe + h - c -)> (16) 

U (1) = f(fioE' "UwW* + hx -) ; (I 7 ) 

u {0) = ^J2' & l + A-^>*'> ( 18 ) 
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where the prime over the summation is used to indicate the omission of the zero mode for the summation. In the 
large N approximation, we keep the term only. The rest contributes only in higher order corrections. 

Next we eliminate the zero mode creation and annihilation operators as follows [38]. Define a self-adjoint operator 
A) by 

00 = (N + lT 1/2 a , N = aja . (19) 

This operator /3 and its hermite conjugate satisfy ftofto = 1 an d A) A) = 1 — fto, where flo = | A^o = 0)(Nq = 0| is the 
projection operator onto the state No = 0. In the presence of the homogeneous BEC for the zero mode, we can safely 
exclude the state No = 0. We then approximate /§o/3o — 1, i-e. (Jo is an almost unitary operator, and 0o,Po\ — 
holds. More precisely this approximation holds between matrix elements when we compute expectation values. Then 
Oq = V~NoPo an d a = /3 \/No allow us to eliminate the bare zero mode operators a J and a . We introduce a new 
set of the creation and the annihilation operators and a k for non-zero momentum by 

4 = OfcA). = Pl&k for k ^ 0, (20) 
which satisfy the equal time canonical commutation relations : 

[dfe (*),<%(*)] = 6 k ,k' - a k >alfl ~ 6 k ,k', ^ 
[&k(t),a hl (t)] - [a{(t),&l(t)}=0. 

The vacuum state of the Fock space Ho for the new set of the creation and the annihilation operators is the same as 
before, and it is defined by 

dfe|0) = for all k ^ 0. (22) 
It should be noted that these composite creation and annihilation operators have the following properties : 

&i\N ,N lt --- ,N k ,---) = y/N^\No-l,N lt --- ,N k + l,---}, (23) 

and 

a k \N ,N lr -- ,N k ,---) = v^|JVo + l,JVi,--- , JV fc - 1,- ■ ■), (24) 

where N k is the number of particles occupying the mode k in the plane wave basis. In other words, the total number 
of particles is conserved after we apply these operators. 

We express the Hamiltonian (7) in terms of these new operators, N , and N, 



Ha = W (N - iV )+ 2^^ + fN )ala k + ^No + 2^No + l}^ ^d T _ fe + h.c). (25) 
Similarly for the interaction Hamiltonian (8), 



Hi = -No + -(\/No + 1 J2' Pk(t)&i + h.c.) +yJ2' ~P« (*)4+*'«fc- ( 26 ) 

We note that the total number operator N = N + a\a k commutes with both Hamiltonians (25) and (26) at 
this point. Therefore we can choose the basis in which TV is diagonalized with given value N. This gives the condition 
No + a\,a k = N which allows us to eliminate the zero mode occupation number operator No- Using an expansion 
for the square root function and keeping up to the second order in the new set of creation and annihilation operators, 
we obtain our approximated Hamiltonian as 

H tot =Eo+ ]T'(e fe + gn)a{a k + ^ ^'(4 At - fe + h.c.) + + ^.c). (27) 

Here E = gnN/2 — gn/2 + An is a constant term, and we neglected terms of order of iV _1 . 
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B. Diagonalization 

Let us implement the Bogoliubov transformation to diagonalizc the bilinear terms in the Hamiltonian (27). The 
required canonical transformation U = exp(iGB) is generated by 

GB = ^e k &i&l h + h.c, (28) 

with 

6 k = tanh"^- ^— ), frw k = y/e k (e k + 2gn). (29) 

Bogoliubov's excitation (the "bogolon" ) is created and annihilated by the operators 

b{ = e lGs a{e' lGB = a{ cosh k + d_ fc sinh 9 k , (30) 
b k =e lGB a k e' lGB = a k cosh 6 k + a ] _ k sinh 6 k . (31) 

These operators satisfy the same canonical commutation relations as before. The new vacuum state of the Fock space 
Hb for the bogolon is defined by 

b k (t)\0 B )=0 for all k ^ 0. (32) 

This new vacuum contains infinitely many modes created by the original creation operators a k 

|0 B ) =cxp(zG B )|0) (33) 
= cxp{ Y^l- ln(cosh0 fe ) + a£at fe tanh0 fc )]}|O). 

It is well known that two Fock spaces H.q and H. b are unitarily inequivalent to each other because they possess infinite 
degrees of freedom. The above formal expression (33) loses its meaning in the thermodynamic limit. Therefore we 
need to specify the vacuum state corresponding to each Fock space [43, 44]. To see this point more clearly, we evaluate 
the inner product of two vacua |0) and \0 B ): 

(0|0 B ) = exp[-5^'ln(cosh0 fc )] - cxp[-V J -0^ ln(cosh^)]. (34) 

Obviously this results in (0|0_b) = for the V — > oo limit, meaning |0_b) cannot be written as a superposition of 
vectors belonging to the Hilbert space Ho- 

After the Bogoliubov transformation generated by (28), the total Hamiltonian (27) becomes 

H B = E' + ]T' hu k bl b k + J2' Uk{t)b{ + h.c.). (35) 

E, 
(hu! k — e k — gn)/2 is the ground state energy without impurities, and 



M ] = " x \lj^tr k p^- (36) 



In the thermodynamic limit, the second term in E can be evaluated as follows. 

1 v-'^ x V f d 3 k 

- (fi^k ~e k - gn) = — J -—3 V™k - <* - gn) (37) 
f d 3 k 1 64gN(na s ) 3 / 2 

-- N —frJ + — 150? ■ (38) 

Obviously, the first term diverges in eq. (38). This divergence can be understood as an artifact of the approximation for 
the general two-body interaction V(x — x') by the delta function g6(x — x'). As a consequence, the Fourier transform 
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of interaction is constant g for all k. Therefore, all wave numbers including the ultraviolet region contribute with the 
same weight. The standard way to deal with this divergence is to renormalize the coupling constant g by 

^h 2 a s Mg 2 f d 3 k 1_ 

M ~ 9 h J (2tt)3 fc 2 ' ^> 

This renormalization is done with the term gnN/2 in Eq, and corresponds to including the second order Born 
approximation for the scattering length a s [45]. 

The first order impurity effects in our model are described by the model equivalent to the generalized time dependent 
van Hove model. We will explain later that the factor l/\/~N due to BEC introduces the significant difference and 
gives rise to the macroscopic degree of freedom. In this sense it is different from the van Hove model. 

The Heisenberg equations for the bogolon operators are 

ihosi = [b{,H B \ = -hio k b{ - fm, (40) 

ihd t b k - [S fc , H B ] = f^kb k + f k (t). (41) 

We solve this differential equation with the boundary condition such that the operators evolve without impurities 
asymptotically at the remote past t = — oo. The solution is 

= + «(*)> (42) 

b k {t) = ti£{t) + 4> k {t). (43) 



Here 4>k{t) is a c- number function : 



,(t) = —e-^ [ dt 1 hit'y^' (44) 



inX Itk , .,, - i0Jkt 



h V Nhuk 



I k {t)e- W *\ (45) 



h(t) - / dt' pkity^ . (46) 



where the integral Ik(t) is defined by 



Alternatively we can directly diagonalizc the Hamiltonian Hb using the time dependent unitary transformation 
Uc(t) = exp (iGc(t)) a s will be shown below [46]. 

In the Heisenberg picture the time evolution of operators 0{t) is governed by the Heisenberg equation ihO(t) = 
[0(t), H{t)\. (Explicit time dependence is not assumed for 0{t).) The unitary transformation U{t) = exp(iG(t)) 
for the operator 0(t) introduces a new operator 0(t) = UOW . Time evolution of the new transformed operator 

0{t) is then given by the Heisenberg equation ihO(t) = [0(t), H(t)]. The transformed Hamiltonian is H{t) = 
U(H — ihWU)U\ where the left hand side is expressed in terms of new transformed operators. In our case we 
implement the time dependent unitary transformation Uc{t) generated by 

G c (t) = -i^ cj> k (t)b{+h.c. (47) 
k 

4>k{t) is a c- number time dependent function which will be determined later. This transformation shifts V k (t) and 
k{t): 

KHt) - e^'bi (t)e- iG ° = b{(t) <j>* k (t), 
&fc(t) = e lGc b k (t)e~ lGc = b k (t) - fc (t). 

To obtain the transformed Hamiltonian in terms of b^(t) and &^(f), we need to evaluate K(t) = UqUc- The formula 

ex P HG)-|exp(iG) = iG + i 2 {G,G}/2\ + i 3 [[G,G},G}/3! + ■ ■ ■ , (49) 
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after simple algebra yields, 

K(t) = {iu) k cf) k + <pk)b k + (iuk<t>* k - <i>k)bk + ^[{i^k<t>k - 4>k)4>k + (i^fe^fe + 4>k)<t>V[- ( 50 ) 
The transformed Hamiltonian Hsit) is 

H B (t) = U c (H B (t) - ihK{t))ul (51) 

= E (t) + hu k b?{t)i% (t) = H in (t), (52) 

provided that the c-number function 4>k{t) satisfies a simple differential equation: 

ih<fi k (t) = hu k 4> k (t) + f k (t). (53) 

The solution with the initial condition <fi k (t = — oo) = is given by (44). In the diagonalizcd Hamiltonian we defined 
the new ground state energy Eo(t) by 

E (t)=E' +Y^Re(f*k(t)Mt))- (54) 

The asymptotic creation and annihilation operators denoted by b k (t) and b k (t) respectively satisfy the same 
canonical commutation relations, and time evolution is governed by the Heisenberg equation with the diagonalizcd 
Hamiltonian H in , i.e. b%\t) = e^b^ and b™{t) = e - iuJkt b%. Hence, 



The energy spectrum huj k is that of the gapless excitations characterized by ivk — kc for a small k, where c = y/gn/M 
is the speed of sound. Fig. 1 shows the characteristics of Bogoliubov's spectrum. We remark that the spectrum ujk 
and the speed of sound c are the same as in the original Bogoliubov model without impurities. Therefore the motion 
of impurities does not affect either w k nor c in our approximation. Another remark regarding Bogoliubov's spectrum 
is that the general two-body interaction V(\x — x'\) instead of the contact interaction gS(x — x') produces different 
curve as shown in Fig. 2. This Fig. 2 is plotted for the finite height potential 

V{\x-x>\) = { V °() X } <ro } (56) 
1 U \0 (|*| >r ), V ' 

with Vo and r constants. The contact interaction approximation is understood as the limit of this potential. In the 
general two-body interaction case, the excitation spectrum is 



hu k = ^e k {e k + 2nV(k)), (57) 
where V(k) is the Fourier transform of the potential 

V{k) = J d 3 x V(|x|)e- ifc -". (58) 
For the above particular choice of interaction (56), 



. i 2Tmh 2 r V f sinkro 



*** = fl+ M \~kW ~ C0S fcro ) ■ ^ 

This excitation spectrum has also phonon-like behavior for small k, i.e. ftw k ~ Kkd with d = \/2irnVorf ) /(3M). 
Interestingly by taking into account the finite distance interaction, Bogoliubov's excitation contains the roton-like 
excitation spectrum as seen in Fig. 2 [47]. 

The vacuum state of the Fock space Hi n for ti k n ^ and b 1 ^ is defined by 

6j c n |0 in )=0 for all fc t^O, (60) 
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FIG. 1: Bogoliubov's excitation spectrum hu>k = sj ek(ek + 2pn). hu>k/(Mc 2 ) is plotted as a function of hk/(y/2Mc). Two 
asymptotic curves are also plotted. The dashed line represents hkc (fc£ <C 1). The dashed curve represents e k + Mc 2 (fc£ >• 1). 




i i a n ^ 



FIG. 2: Bogoliubov's excitation spectrum ftcjfc = yetfe* + 2nV"(fc)) for (59). hu} k /e kQ is plotted as a function of k/ko with 
ko = 1/vq. 

and this vacuum represents the homogeneous condensate. The relation between two vacua of Fock spaces He and 
Hi n is formally expressed as 

|0i») =exp(iGc(t))|0 B ) (61) 
= exp[^' ( - l|0 fc (t)| 2 + M0hWb). (62) 

The homogeneous condensate vacuum |0; n ) is time independent since the time dependence of the annihilation operator 
bfe(t) is simply exp(— iukt). Hence the Hilbert space of bogolons is a one parameter family of H[ n characterized 

by the time t. 

The physical interpretation of the above result is quite simple. The motion of impurities creates the time dependent 
coherent states. The dressed Bogoliubov's excitations do not annihilate the homogeneous condensate vacuum: 

&fe(*)|O in ) = fc (t)|0 in >. (63) 

III. OCCUPATION NUMBER AND ENERGY DISSIPATION 

One of our main effects described by our model is the effect of energy dissipation due to the motion of impurities 
in the homogeneous BEC. From our model we find that even though the systems is in the ground state where no 
elementary excitation is present, there is energy transferred from the motion of impurities to massive bosons. In 
other words, the motion of impurities will dress bogolons in such a way that bogolons will not see the homogeneous 
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condensate as a vacuum. Since we do not take into account any back reaction on the impurities, this transferred 
energy is the energy in order for the impurities to move along given trajectories. It means that we need to feed this 
amount of energy to keep the impurities moving. More realistic model for the dissipation effects in BEC should be 
done by considering quantum impurities and taking into account the dynamics of them. 

We first evaluate the occupation number for bogolons with respect to the homogeneous condensate. This number 
counts the emitted bogolons accompanying with the motion of impurities in BEC, 

h k (t) ee (o in |St ( t )S fe (t)!o in ) = \Mt)\ 2 = ^|IM*)I 2 - (64) 

Fluctuation in the occupation number Ah k (t) is 

An fc (t) = l(0 in \(b{(t)b k (t)) 2 \0 in ) - ((0 in \bl(t)b k (t)\0 iD )) 2 }^ 2 (65) 
= \Mt)\- (66) 

We are interested in evaluating the occupation number for the t — > oo limit, which covers the whole trajectories of 
impurity motion. 

We next evaluate the expectation value for the number operator N k {t) = d\{t)d k {t) = a\{t)a k {t) (k ^ 0) with 
respect to the homogeneous condensates |0j n ) as 

<o in |^)|o in ) = IC-^f - 1) + ^ + Ulfc(*)l a + - DI*-*MI 2 + £-J€(t) - 0- fe (t)| 2 . (67) 

The first term is the well known term in the original Bogoliubov argument, and the rest are contributions due to the 
impurities. The depletion of the condensate d(t) due to the quantum fluctuation can be evaluated by summing over 
all modes and divided by the total particle number N: 



*» - iV ^ + s + &im - ♦-.WP). <«> 

Since \<j> k \ 2 has an additional factor l/N, the motion of impurities does not contribute to the depletion of the homoge- 
neous condensate in the thermodynamic limit. Therefore, in the large N limit, the homogeneous condensate is stable 
against a small number of external impurities, i.e. when the number of impurities is much less than the total number 
of bosons. In contrast, real experiments have always the finite number of particles and a finite size. Therefore, there 
are effects on the depletion of condensates due to the motion of the impurities. 

The dissipated energy £ k {t) for a given mode k is obtained by multiplying n k (t) by the excitation energy hu k} i.e. 
£k(t) = fkOkn k (t). Therefore the total dissipated energy £{t) is given by summing over all modes, which is expressed 
in the thermodynamic limit as 

E{t) = -> J J^ ek \I k (t)f. (69) 

In particular for the t — > oo limit, taking into account the complete trajectory of impurity motion, we obtain rather 
a simple expression, 

Stat - lim £(t) (70) 

t— *oo 



n\ 2 f d A k 



dui 5(iv - L) k )k 2 \p c (k,u))\ 2 . (71) 



2M J (2tt) 3 

Here p c (fe,w) is the Fourier transform of impurity distribution p c (x,t) with respect to space and time coordinates : 

p c (k,u)= I dt ( d 3 x p c (x,t)e- lk - x+wt . (72) 



IV. EMERGENT MACROSCOPIC FIELDS 



The coherent states play an essential role in representing classical fields from quantum systems in various areas of 
physics. In quantum electrodynamics for example, the classical Maxwell equation can be derived as the Ehrenfest 



11 



equation with respect to the coherent state. In our model one of the main results eq. (63) is best understood form 
this point of view. To this end we express operators in terms of incoming field creation and annihilation operators 
b 1 ^ and S£. 

We first express the field operator ip(x,t) in terms of incoming fields. For this purpose we need to know the 
time dependence of operator $o(t) defined by eq. (19). Using the Hamiltonian (25, 26), it is easy to show that the 
Heisenberg equation for this operator is ihdt$o{t) = gn$o(t) within our approximations, and hence 

$ (t) = e -i9nt/h pin = g-iMc 2 */^. (73) 

Using the Bogoliubov transformation (31) and the solution (43), 

${x,t) = ^e- mc2t ' h [^+ ^^'(cosh# fe b k (t)e ik - x -sinh0 fe S£(t)e- ifc -")] (74) 



V 

= ^e- iMc2t / h (M^t)+Mx,t)). (75) 
Here the c- number field tp c (x,t) and q- number field ip q (x,t) are defined by 

VcOM) = Ml - ^ E' -^— {haJkIk{t)etk " kt *kl* k (t)e- ik - x+i " k % (76) 
j> q (x,t) = ^E' + £k + hLUkb k elk " kt - Van + e k - h7 k ^e- ik - x +iw k t), (77) 

with Ik{t) eq. (46). A similar expression for ^(x 7 t) can be obtained by taking the hermite conjugate of the above 
result. 

We next evaluate the density operator p(x, t) and the current density operator j(x, t) defined by 

p{x,t) = ${x,t)i){x,t), (78) 
iOM) - ^[^(x.tXV^x.t)) - (V^(x,t)Mx,t)]. (79) 



(80) 



We express these operators in terms of incoming fields in the same manner as before, 

«-.«) -4r M mt)er ""' + «■> + 7v £ M &<,)e ~""~ + h • c • ) • 

i(*.*) = 7 E wV^ ( ^ (<)e +c - c - ) + w E 2mV^t (6;: We +h - c - } - (81) 

These equations satisfy the continuity equation : 

ftp(x,t) + V-j(aj,t) = 0, (82) 

as required by the self-consistency of our approximations. In the thermodynamic limit, these operators have non- 
vanishing expectation values with respect to the homogeneous condensate |0; n ). We now define the induced density 
due to the motion of an impurity by 

p{x,t) = (0 in |p(;M)|0in) - (0in|/3(x,i)|0 in )|A = 0. (83) 

In the homogeneous BEC case, this induced density accounts for the deviation from the constant density n. The same 
expression for the induced current density is defined. Then eqs. (80,81) yields 

/rl 3 k k 2 
J^^ k {lIk{t)e ^ k ' X+lUJkt + C - C ^ (84) 

/rfik k 
j^2 {im e- ik - x+i ">< t + c.c.), (85) 

where we denote the ratio of two coupling constants g and A by k = X/g. These induced density and current density 
also satisfy the continuity equation : 

d t p(x,t) + V-j(x,t) = 0. (86) 
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Although the Bogoliubov excitation frequency uj k contains h, the cancelation of h in the above expressions should be 
noticed. 

It is also convenient to obtain the Fourier transform expression in terms of p c {k, u) (72). Using the formula 



lim / dte twt = dt e^-^ = , (87) 

r^°o,/ J a w + ie 



/' 

T^oo 

with e a positive small number which goes to zero at the end of calculations, p(x,t) and j(x,t) become 

Fourier transforms of p(x,t) and j(x,t) are written using the Fourier transform of the retarded Green function, 
D ret (w, uj k ) = 1/(uj 2 -cj 2 + isgn(w)e), 

p(fc, w) = Kc 2 k 2 p c (k, w)D ret (w, Wfc), (90) 
j(fe,w) = K,c 2 kwp c (k, w)D ret (w,Wfe). (91) 

The induced density satisfies the following differential equation: 

(<9 2 - c 2 A + c 2 ^ 2 A 2 )p(a;, i) = kc 2 Ap c (s, i), (92) 

where A = V 2 is the Laplacian and £ = h/(2Mc) is the coherence length of the homogeneous condensate. We can 
also obtain the differential equation for the induced current j(x,t) in the same manner. Assuming the distribution 
of impurities p c (x,t) satisfies the continuity equation d t p c (x,t) + V • j c (x,t) — 0, e.g. for a point-like impurity 
p c = 5(x — £(t)) and j c = £(t)5(x — C(*)) 5 we have 

[d 2 - c 2 A + c 2 ^ 2 A 2 )j(x, t) = -nc 2 Vd tPc {x, t) (93) 

= kc 2 V(V -J c (x,t)). (94) 

Lastly we evaluate the induced stress energy momentum complex T^tJx, t) (p, v = 0,1, 2, 3) defined by 

f^ } (x,t) = (0 in |T"J„ ) (x,t)|0 in ) - (0 in |f^ ) ( a; ,i)|0 in )| A= o. (95) 
The stress energy momentum complex T^(x,t) of the Bogoliubov model is obtained from the Lagrangian density 

c = y (ftw - d t 4>H) - ^diftdd - ffi&w, (96) 

where the space-time coordinate dependence (x,t) is for simplicity omitted for ft(x, t) and ip(x,t). When the 
action J d 3 x J dtC is invariant under an infinitesimal space-time translation, there exists the stress energy momentum 
complex of four conserved Noether currents [53, 54]. These four conserved currents are defined by 

Fir ~ ~ Ftr 

where <9 — d t and summation of repeated indices is understood. For the Bogoliubov model, 

f\ 0) = -^(8^84 + 8^8^), 

T° (i) = y(^0iV- ( 98 ) 
h 2 



fi U)=-2M^ d ^ + d ^ d ^ 



+ Sij[~ (ft 84 - d^H) + ^d k ft8 k 4> + 9 -ft 
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These four Noether currents satisfy the continuity equation c^T^ n =0 [y = 0, 1, 2, 3), which results in conservation 
of the energy {v — 0) and the momentum {v = 1,2,3). For the general two-body interaction case, the stress energy 
momentum complex does not conserve locally, but it does conserve globally. 

Substituting eq. (74) and its hermite conjugate into eqs. (98), it is easy to show that the Fourier transform of the 
induced stress energy momentum complex, 



T"i v) (k,u) = J d 3 x J dtT" {v) (x,t)e 



(99) 



is 

f° {0) {k,w) = KMc 4 k 2 p c (k,oj)T> iet (w,co k ), (100) 

T l ( Q) (k,uj) = nMc^kiOjpdk, u>)~D let (u),(j k ), (101) 

T° (i) (fc,w) = -KMc 2 /ciwp c (fc,w)D rct (w,w fe ), (102) 

T^fc.o/) = S ljK Mc 2 LU 2 p c (k,uj)-D Ict (Lj,Lj k ). (103) 

From the above result, we see that it is natural to define the effective Minkowski metric tensor 

V, = diag(c 2 , -1,-1,-1), (104) 

with c the speed of sound, and its inverse 

rf v = diag(l/c 2 , -1,-1,-1). (105) 
Then we raise and lower the indices to get the symmetric stress energy momentum complex defined by T^{k,u) = 

f oa (k,iu) = KMc 2 k 2 p c (k,uj)D rct (uj,Lj k ), (106) 

f M {k,uo) = KMc 2 fc^p c (fc, W )D rct (w,^), (107) 

T 0i {k,uj) = KMc 2 ki0jp c (k,u)D ret (u,uj k ), (108) 

T ll {k,uj) = 5 i0 KMc 2 uj 2 p c {k,uj)T), c t{ijJ,Lo k ). (109) 

This induced stress energy momentum complex is conserved locally in our approximation 

d lx f llv {x, t) = (y = ~ 3). (110) 

Several remarks are in order. We first observe that for very large length scale, much greater than £, we can 
approximate eq. (92) as 

V^d^d u p(x, t) = kA Pc (x, t). (Ill) 

Thus bogolons propagate following the ordinary sound wave equation in barotropic fluids with the speed of sound 
c regardless the details of the motion of impurities. We emphasize the non-triviality of this result. Although the 
homogenous BEC maintains its condensed state, where the depletion of condensate is very small in the large N limit, 
there arises macroscopic change in its density (0i n \p(x, i)|0; n ) = n + p(x, t). Especially, any small disturbance in the 
homogeneous BEC at some space-time point will produce the propagation of bogolons inside the space-time light 
cone as represented by the wave equation (111). In this large length scale, the left hand side of (111) respects the 
Lorentz symmetry rather than the original Galilean symmetry. This effective symmetry becomes more exact as we 
go to a larger length scale. However, this classical behavior will not hold when we reach the high energy scale or the 
very small length scale £, where we cannot neglect the terms depending on the coherence length £. This situation is 
familiar: quantum hydrodynamics breaks down at certain length scales. There, the impurities start to see the atoms 
rather than phonon-like low energy excitations. 

From the induced stress energy momentum complex, we have f°, , = Mc 2 p, i.e. the deviation of the energy 

density of this fluid is equal to Mc 2 x (the induced density). We also have the relation = T 0t c 2 , in other words, 

(energy density flux) = (momentum density) xc 2 . The famous Einstein relation holds true, with c the speed of 
sound. Lastly, we point out that the induced stress energy momentum complex T fJ,u (x,t) is completely symmetric, 
i.e. T^ v (x,t) — T ufl (x,t). Here we used the inverse of the effective Minkowski metric tensor rj^ to raise the index: 
T^ v (x : t) = rjVT^tJxjt). This is a non-trivial result. This suggests that relativity could indeed emerge as the 
consequence of the existence of an underlying constituent. 
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V. EXAMPLES 

We will consider several examples to illustrate our model and its general solution. In order to discuss the analogy 
to special relativity, we use the relativistic form of trajectories. (In order to discuss realistic experimental settings, we 
need to use the non-relativistic trajectories instead. Moreover, we should use finite times, volumes, and trajectories 
conforming to the experimental situations.) The speed of sound plays the role of the speed of light. We are interested 
in the occupation number (64) and the total dissipated energy (69) for the t — > oo limit. Therefore our main concern 
is to evaluate the integral (46) for p c (x, t) = 8(x — C(i)), with a trajectory of point-like impurity. 

J fc (oo) = J dt J d?x p c {x,t)e- lk - x+tu > kt (112) 
dt e *"**-ife-C(*). (113) 



= / 

j — ( 



A. Method of Stationary Phase and Landau's Criterion 

Since it is impossible to carry out the above integral (113) in general, let us first examine it using the method of 
stationary phase. 

We are interested in the behavior of this integral for high momenta where the condition ojk S> fee holds. Then the 
integral oscillates rapidly, so that we can approximate it using the method of stationary phase. Rewrite the integral 

as 

J fc (oo) = / dte™^, (114) 



where g{t) = t — k ■ C,{t)/u)k = t — k ■ £(f)/(c^/l + (fc£) 2 )> w i tn & = k/fc a unit vector in the direction of k. The 
stationary phase point is obtained by 



g'(r) = 0^k- C(r)/c = y/T+ (H) 2 (> 1), (115) 

and the integral is approximated as 

Jfc(oo) ~ ^e^J-tt-Cfa) f dt e -*fc-C(Tj)(*-Ti) 2 /2 (H6) 



_ ^^Tj-ik-cw l—=Z . (117) 

j y ife ■ C{ T j) 

Here the summation should be taken over the stationary phase points obtained by (115), and k ■ C(tj) ^ is assumed. 
Otherwise the integral is approximated as zero. In order to give rise to the manifest energy dissipations, the argument 
of the exponent (117) has to be imaginary. However as long as the velocity of the impurity is less than the speed 
of sound c, no real r exists satisfying (115). This means any subluminal motion below c will produce exponential 
suppression factors resulting in much smaller dissipation effects. (In addition we need a condition lm(g(rj)) > 0, 
otherwise the integral grows exponentially.) This general statement is equivalent to Landau's criterion [48]. Indeed, 
if we start from the general two-body interactions V(x — x') rather than the contact interaction as mentioned before, 
the above stationary phase point condition becomes 

fe-C(t)=w(fc). (118) 



2tt 



Here u{k) = \J efe(efe + 2nV(k))/h is the dispersion relation for the general two-body interaction. Thus we obtain the 
critical velocity v c as 

»c - nrin ( ^) ■ (119) 



fe V l fc l 

The stationary phase points are always complex for \£(t)\ < v c . The corresponding statement is that there is always 
an exponential suppression when the speed of the impurity is smaller than the critical velocity \C{t)\ < v c . Therefore, 
we derive Landau's criterion from the first principle based on our simple model. We will give a more quantitative 
argument by considering the motion of an impurity with a constant velocity later. As is clear from our discussion, an 
impurity moving slightly slower than the critical velocity can give rise to energy dissipation. This cannot be explained 
from Landau's argument based on the kinematical consideration. 
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B. Two Stationary Impurities 



As a special case, we first consider two point-like impurities sitting in a BEC. For simplicity we assume the impurities 
have the same mass, i.e. the coupling constant A is the same. The distribution is p c (x) — (6(x — X\) + S(x — x 2 ))/2 
where the factor 1/2 is a convention to make / d 3 x p c (x,t) = 1, and xi, x 2 are constant vectors. It is clear that this 
special case does not give rise to energy dissipation in the t — > oo limit, since Jfc(oo) oc 5(iVk), and u) k = if and only 
if k = 0. However, the ground state energy 

Eo(t) =E' + ^Be(fZ(t)Mt)), (120) 

depends on the distance between the two impurities R = \x\ — x 2 \. 

The limit t — > oo should be taken carefully in this case, since the Fourier transform of p c (x) does not depend on 
time, i.e. f k is time independent. 

V NH 2 LU k T-,ooJ_ T ) 

-n 2 \ 2 e k \p k \ 2 l-cosw fe (t + T) 
= ^2 llm ( 122 ) 

Nh Z LOk T— oo UJk 

_ -n^e k \p k \> 



NH 2 uj 2 ' 

where we use lim-r— >oo(l — cos wT)/w = V(l/w) with V the principal value distribution, and V(l/w) = l/w for i»/0. 
The above limit holds at any time t except t = — oo. Hence the ground state energy is a constant of time, and in the 
thermodynamic limit, 

. MnX 2 f d 3 k 1 gn 2 X 2 f d 3 k 1 n\ 2 f d 3 k cos(k-R) 



h J (2tt) 3 k 2 h 2 J (2tt) 3 lu 2 2 J (2tt) 3 e k + 2gn 

= E ^^2^-^W R • (126) 

The divergent integral, the second term in the second line, should be considered as a term arising from the fact that 
we approximate the interaction Hamiltonian (4) by taking only contact interaction with the s-wave scattering length 
b 8 . We have already encountered the same divergent integral when we calculated the ground state energy E' without 
impurities by renormalizing the coupling constant g. The same procedure is applied for the coupling constant A which 
appears in E' . 

The third term in eq. (126) has the same form as the Yukawa potential and gives rise to an attractive force between 
two stationary impurities. In real experimental situations however, the distance between two impurities is much larger 
than the coherence length, then attractive force decays with the exponential factor. 

On the other hand, if the mass M of boson happens to be very small such that the coherence length has the same 
order as the size of the system L = (V) 1 / 3 , then the effective potential V e s for the impurities can be approximated as 

VcS * R- (127) 

Thus the homogeneous condensate produces a long-range attractive force between impurities. This is the same as the 
Newtonian gravity regardless of details of the impurities. The strength of this attractive force depends only through 
the combination n\ 2 M / (Aith 2 ) . Expressing the coupling constant A in terms of the s-wave scattering length b s , the 
coefficient of (127) is approximated n\ 2 M / [A-kH 2 ) ~ nirblh 2 / M . In this scenario, the impurities with a longer s-wave 
scattering length get attracted stronger. 

We remark that this effective attractive force was also obtained by others by studying a mixture of bosons and 
fermions [22, 50, 51]. There they found that a bosonic background produces an attractive Yukawa type force between 
two-component fermions. 
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C. Constant Velocity 

As the simplest example of the motion of an impurity in the homogenous BEC, let us consider a linear motion of 
one impurity with a constant velocity v. The trajectory is £(t) — vt, then the integral (113) is 

Jfe(oo) = 2tt 8{uj k - k -v). (128) 

To obtain the total dissipated energy in this case, one needs to give a little attention to the limit t — > oo since we 
cannot take a square I k (oo) directly. To this end we first integrate for a finite time interval T, and then we take the 
limit T — > oo at the end of calculations. The dissipated energy becomes proportional to T for large T with the use of 
the formula (1 — cos wT)/w 2 ~ ttT6(w). Hence we obtain the rate of energy dissipation T in the thermodynamic limit 

r= km (129) 

27mA 2 f d 3 k , . 

e k 8(uk - k ■ v). (130) 



h 2 J (2tt) 

Thus T = unless there exists k such that ujk = k ■ v. To proceed further, we rewrite the above integral in spherical 
coordinates and introduce the new variables x = hk/ (2Mc), and y = (3 cos 9, where (3 = v/c. Then the rate of energy 
dissipation reads 



4n\ 2 M 3 c 



r= nhij3 j o dxx 3 j dy 5{y/^+l-y). (131) 
Clearly it follows that (3 > 1 in order to have energy dissipation. Hence the rate of energy dissipation is 

r = ^^-(/3 2 -i) 2 e(/?-i), (132) 

where 0(x) is the step function. Like the Cherenkov radiation, this radiation is emitted inside a cone with an opening 
angle 6 C = cos _1 (l//3), and its momentum cuts off at hk c = 2Mc\J (3 2 — 1 keeping the total dissipated energy finite. 
As was mentioned above, this is the microscopic description of Landau's criterion. We now present the original 
argument by Landau, which is based on purely kinematical considerations. 

Consider any excitation in a fluid, which has atoms of the mass M, has an energy E and a momentum p in a reference 
frame where the fluid is at rest. Then the energy of the excitation E' in the coordinate system where the fluid is 
moving with velocity v with respect to the container, is given by a Galilean boost: E' = E+p-v + Mv 2 /2. Dissipation 
due to the motion of the excitation can happen if the energy in the boosted frame is negative, i.e. E + p ■ v < 0. In 
the optimal case where p is anti-parallel to v, this implies v > E/p. Therefore, for the phonon-like excitation E = pc, 
this condition gives v > c, with c the speed of sound in the fluid. Superfluidity is explained then by the absence 
of excitations for motions with velocities smaller than the critical Landau velocity. While this condition is widely 
accepted, it seems that it has not yet been derived from an underlying microscopic theory as was done in this section. 



D. Constant Acceleration 



The trajectory of a linear accelerated particle along the z-direction is 



CW = ^l + (f) 2 z, (133) 

with a given acceleration a and c the speed of sound. The speed of the impurity is always bounded by that of sound, 
\C(t)\ < c. In the limit c — ► oo, the trajectory becomes quadratic in t as at 2 /2+ (constant). To evaluate the above 
integral (113), we change the integration variable t to a dimcnsionless variable s through t — (c/o) sinhs, and we 
define 

n k = C -^Lu 2 -(ck-z) 2 = ^-^/ S m 2 6 + k 2 e, (134) 

, i ,ck ■ z. , _i , cos# , , . 

a k = tanh- 1 (— — = tanh { -= , (135) 
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where is an angle between k and the z-axis. Then, 

c f°° 

Jfe(oo) = - / rfs coshs exp[iOfe sinh(s — <Tfe)] (136) 

a J—oo 

2c 

= — [it 5(Qk) cosher*, + iKi(Ofe) sinherfc], (137) 

where Ki(a;) is the modified Bessel function of the second kind. The first term in cq. (137) will not contribute to 
the occupation number for the following reasons. The first is that flk = if and only if k = 0. The second reason 
follows from two properties of the Dirac delta distribution; S(f(x)) — J2 i S(x — Xi)/\f'(xi)\, where Xi are the roots of 
f(x) = 0, and x8{x) = 0. The final reason is the presence of the factor et/uJk in cq. (64). Therefore the occupation 
number is 

n k = 4 f," 2 ^ — (Ki(» fc )sinhcr fc ) 2 . (138) 
h A a 2 N Lu k 

The total dissipated energy is expressed as 

Let us examine the above results for the infrared and the ultraviolet region using the asymptotic behavior of the 
modified Bessel function. For this purpose we will take differential forms of the occupation number dhk and the 
dissipated energy d£ to t = hio^dhk in the thermodynamic limit. From eqs. (135,138), 

2nAWW<? Ki(nfc) 2 
dUk = (2nfMa^ k { ^ k ~ ) dM °> (W0) 

with flk eq. (134) and do the element of solid angle. 

Qk becomes very small in the infrared region, and hence the leading term of eq. (140) is 

2™A 2 fccos 2 6> J;J 

dn k ~ ; dkdo. (141) 

(2n) 3 Mc 3 hy/l + /c 2 C 2 (sin 2 6 + k 2 ^ 2 ) 2 

Therefore, the leading term has a singularity around k ~ at the angles — 0, n, i.e. along the direction of trajectory, 
which results in the divergent result for the total number of emitted bogolons. This infrared singularity due to bogolons 
has the same origin as the famous infrared catastrophe in quantum electrodynamics [49]. There, the total number of 
emitted photons due to a motion of accelerated charge particle diverges as k — > 0. However, the physical measured 
quantity, the radiated energy, is still finite. 

To see the comparison between our case and the infrared catastrophe in quantum electrodynamics, we first integrate 
cq. (141) over the angles. Then we find that the leading term of dhk is dk/k, besides a few factors, which diverges 
for k — > 0. However, the total dissipated energy d£tot is finite after we multiply dhk by hwu — hkc. 

We next use the asymptotic form of Ki for large Vl k to get the leading term: 

n\ 2 (2Mc) 3 cos 2 e , hk 2 c^ „ , 
dHk * (2«n*ak* eM-^dkdo. (142) 

Hence the dissipated energy in the ultraviolet region after integrating over the solid angles is 

4nA 2 M 2 c 3 , hk 2 c, „ , „ o , 

^-^^ exp( -W*- (143) 

Thus d£tot asymptotically depends on k through the exponential factor exp[— hk 2 c/(Ma)] [52]. This high momentum 
behavior could be used to measure the number he/ '(Ma), from which we can obtain the mass M of a boson. 

Lastly we notice that two length scales in this example have different orders of magnitudes. One is the coherence 
length £ and the other is l a = c 2 /a. In most real experimental situations £ <C l a , since accelerations available in 
laboratories are very restricted. Then we again use the asymptotic form of Ki to estimate the total dissipated energy 
for the weak acceleration limit, 



nA 2 M 2 c / ha 



--tot 



10ft 3 V ttMc 3 ' 



(144) 



Although the speed of impurities never exceed the speed of sound we expect a small finite amount of dissipation due 
to the uniformly accelerated impurity. 
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E. Circular Motion 



The trajectory of relativistic circular motion on the xy-plane is discussed in Appendix. Thus it is enough to consider 
the following trajectory, given two parameters R and ft > 0. 

( x (t) = R cos fit, C v (i) = R sin fit, ( z (t) = 0, (145) 

where the constant term in ( x (t) is dropped since it does not contribute to the dissipated energy. The integral to 
evaluate is 

/oo 
dt e iuJ kt — ik±Rcos(i~lt — ifr)^ ( 14 g) 

where we express k in cylindrical coordinates as k = (k± cos0, sin0, k z ). A formula 

oo 

e -ixco S e = J- (-i)V^(z), (147) 

£=-oo 

with 3e(x) the Bessel function, expresses the above integral in a series as 

/oo 
dt e ^ kt+imt (148) 
-oo 

= 2tt i~ e e Ulp J e (k ± R) S(u k - tfl). (149) 

e 

It is clear from the last expression that the above integral vanishes unless uu = tfl for some positive integer t. When 
we square (149) we need to take care of the delta distribution properly in the same manner as with the linear motion 
case. The occupation number for a finite but large time interval T has an asymptotic form: 

2 \ 2 00 

h k = 2nT^^-J2 - ^){h{k±R)f. (150) 
Hence the total dissipated energy per unit time is 

r = lim (151) 



_ 27mA 2 ^ r d 3 k 
1=1 J 

To carry out the angle integral we use the formula 



\ 2 00 r j3 i, 

W~ £ J Wf ek 5{LOk ~ m ^^ R )) 2 - ( 152 ) 



(J £ (x)) 2 = -!- / dipe 2Uv 3 Q {2xsinip). (153) 

2lT J-n 



After several straightforward steps, we get the £th mode for the rate of energy dissipation T: 

mo 2 



T * = F ° 1 rf^ E^+i^fi), (154) 



where To = n\ 2 Mfl/ (2nh 2 R) is a constant with dimension of energy per unit time, and kt is 



fc^=y^/l + 2(^) 2 -l (£=1,2,-..). (155) 

The total dissipated energy is given by summing over all modes I. 

In order to estimate the large i behavior of (154), we use the method of stationary phase as discussed before. The 
stationary phase points are given by 

k ± Rfl sm(flTj - (j>) + uj k = 0. (156) 
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Limiting the case RVl/c < 1, i.e. k±Rtt < Uk, the tj are 

-3 i Ldk 



nr 1 -^=(- + 2j)TT±i\cosh-\ I -^)\ (j = 0,±l,---). (157) 

We need to choose Im(rj) < to satisfy the condition that the square of the integral will not diverge for k — > oo. 
Then, following the procedure discussed before and deforming the contour properly, we get the asymptotic expression 
of the integral Ifc(oo) : 

Ifc(oo) ~ exp[iwfeTj — ifcj_i?cos(f2rj — <f>)] J dt exp[— Q\J — (kj_RQ) 2 (t — Tj) 2 ] (158) 



j = -oo 



1/2 ; 1 r. . I „„„u-l/ W fe 



= ( - - ,, = =) 1 /^exp{--K|cosh^( 



3 
2 



9 00 

+ ^2-(fc ± i?O)2- iWfe (0+^)]} ]T e 2 ™ /0 (159) 



j = -oo 



y^-(*iffl) 2 + - (fc±^) 2 

x ex P {-I[^2 - (fc ± i?fi) 2 - iu> k (4> + |tt)]} 5(w fc - £Q), (160) 

where we used formulae ^2JL_ 00 ex Pfc#) = 27i"X]fc-oo <5(^ — 2^7r) and cosh -1 x = ln(x + \/a; 2 — 1)- Thus wc have also 
obtained the same delta distribution dependence as eq. (149). Taking a rough approximation \/u'f, — (kj_RQ) 2 ~ Wfc, 

7 fe(°°) - E ^y ^^^ - «»■ ( 161 ) 

With this result and the relation I ~ hk 2 /(2MQ) for large f , we find an exponential dependence of the rate of energy 
dissipation for large k as exp[— 7i/c|/(Mf2)] . This shows similar behavior as the constant acceleration case where 
dependence was found as exp[— hk 2 c/{Ma)]. In most experimental settings the linear acceleration is so weak that 
c/a I/O. And hence we expect in the circular motion case, it is easier to estimate the mass M of bosons by 
measuring the dissipated energy for high momenta Hk. 

We remark that circular motion in BECs may be easier to realize in laboratories compared to the constant accel- 
eration case. One reason is that we have two parameters R and f2, which when combined can produce observable 
dissipation effects. The other reason is that circular motion can be realized with additional trapping potentials which 
are necessary in current BEC experiments. In this context it is also interesting to study the induced density and 
current profiles (88,89) to examine the validity of our model. 



VI. DYNAMICS OF IMPURITIES 

We will consider back reactions on the impurities in this section. To this end we regard the impurity as a single 
quantum mechanical particle with the canonical momentum p(t) and the position q(t). They satisfy the canonical 
commutation : 

[quPj] =ihS i: j = 1,2,3). (162) 

We add the Hamiltonian for the impurity H imp = p 2 /2m with m the mass of the impurity. Then the total Hamiltonian 
of the system is 

2vt2 



H(t) = J d 3 x ^( X ,t)(-^)i,(x,t) 



+ | / d 3 x ip^( X ,t)^(x,t)iJ)(x,t)i]j(x,t) 



+ \j d 3 x S(x - q(t))ft(x, t)i>(x, t) + . (163) 
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Here we assume that the impurity is moving without an external potential. The dynamics of the impurity with 
additional potentials can be found by adding suitable a potential term for the impurity in the above Hamiltonian. 
Following exactly the same procedures as before, we obtain the Heisenberg equations of motion : 

ihd t b k (t) = hw k b k (t) + fk(t), (164) 

-ihd t b{{t) = fiw k i>l(t) + fl{t), (165) 
ihd t q(t) = ihp(t)/m, (166) 

ihd t p(t) = ' Hk (bk(t)ft(t) - h.c), (167) 
where f k (t) now depends on the position operator q(t), i.e., 



Nhco k 

Thus we can obtain a closed equation for the position operator q(t) as 

n\ 2 1 x - hke k r 



f k (t) = n\J-^— exp(-zfe • q(t)). (168) 



m 9(*) = -^-^X] ^7 exp(ig k (q(t);t)) dt' cxp(-ig k {q(t');t')) + c.c. 



(169) 



where gk{q(t);t) = oj k t — k -q(t) and the boundary condition is chosen as t = t . To derive this equation, we take the 
expectation value with respect to the homogeneous condensate, and we drop the hat for the position operators. 

Although it is unlikely that we can obtain the analytical solution to the above equation (169), this equation will be 
useful to study the dynamics of the impurities numerically. 

VII. ALTERNATIVE MODEL 

So far we have only considered the local density-density interaction between the impurities and massive bosons. 
In this section we provide an alternative interaction model which is Galilean invariant. We first review Galilean 
symmetry in the non-relativistic quantum field theory [54, 55] . 

A. Galilean symmetry 

A Galilean transformation with a given velocity v is defined by 

t^t' = t, (170) 
x — > x' = x - vt. (171) 

Under this Galilean transformation, the partial derivatives d t and <9j are transformed as 

dt^d'^dt + Vidi, (172) 
d[ -> d t . (173) 

The quantized field t(^(x,t) is transformed as 

$(ss, t) -» $'(x', t') = e iS ^ x ^/ h i,{x, t), (174) 

where S v (x,t) = -Mv x + Mv 2 t/2. 

The Lagrangian density C for the Bogoliubov model, 

C = ^(ftdtj, - dtftj) - ^diftdd - (175) 

is invariant (scalar) under the Galilean transformation. Then the corresponding conserved current N^Jfi — 0, 1, 2, 3 
and i = 1,2, 3) is given as 

N° (l) (x, t) = x l J°(x, t) + t f° (t) (x, t), (176) 
N\ j) {x,t) = x i P{x,t)+tf\ j) {x,t). (177) 
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Here J°(x,t) — Mp(x,t) is the mass density, J l (x,t) — Mf(x,t) is the mass current density, and the stress energy 
momentum complex T^(a;, t) is defined by eqs. (98). For convenience we list them below : 



J°(x,t)=Mft%l>, (178) 

I 1 



J i {x,t) = ^ftd i $-d i ft$), (179) 



f ° (0) (as, t) - ^diftdrf + § (180) 

^(O)^.*) = -^(d t ftW + &iftW), (181) 
f° (i) (a ! ,t) = ^(^^-^% (182) 

+ %hy W^V* - ^V> f V0 + ^M 9 ^d k ^p + (183) 
It is straightforward to show that the current satisfies the continuity equation : 

d t N° {i) + d j N\ i) =Q{i = l, 2, 3). (184) 

Therefore, the generator G 1 defined by 

G* = 1 d 3 ziV (i) (x,i), (185) 

is time independent, i.e. dG % jdt = 0. The commutator between the field ij}(x,t) and the generator G l is 

{4>(x, t),G l ] = (Mi 1 + 1 ihdi)i){x, t). (186) 

B. Current-current interaction 

We next look for a possible interaction with the impurities based on the Galilean symmetry For the purpose of 
symmetry considerations, we first consider interactions between massive bosons and a quantized field ^(x,^ for the 
impurity. The Galilean transformation is also applied for the impurity field *$>(x,t). The basic idea is to add possible 
Galilean invariant terms to the Lagrangian density. The density p = ^ip is obviously a scalar under a Galilean 
transformation. Hence, the simplest interaction is 

Ci = -Xppimp, (187) 

where pi mp = is the density of the impurity field. Taking the classical limit, where the density of the impurity 
is replaced with the Dirac delta distribution S(x — C(*))j we recover the model proposed in Sec. II. 

The next possible interaction is naively expected as the current-current interaction, i.e. C'j = —X'j-ji mp . However, 
this interaction is not allowed within Galilean invariant theory. We note that the current j is transformed under a 
Galilean transformation : 

j(x, t) -» j'(x', t') = j(x, t) - vp(x, t). (188) 

Therefore, the term j - jimp is not a scalar under a Galilean transformation. In order to form a Galilean scalar, wc 
look for a quantity j±{x, t) which transforms as 

j 4 (x, t) -» j' 4 (x', t') = j 4 (x, t)-v -j(x, t) + ^v 2 u(x, t). (189) 
Then, it is easy to show that the following combination 

3 'Jimp PJiirnp j4Pimp: (190) 
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is invariant under a Galilean transformation. We make the following ansatz for Ja(x, t) : 

k(x, t) = - A 5 (^t( aj) t ) V 2 ^(x, t) + V 2 ^(x, t)i>(x, t)). (191) 

Another possible form is J4 = h 2 do\) / (2M 2 ) . The physical meaning of j4(x, t) is clear, if we multiply by the mass 
M. Mj4(x,t) is the kinetic energy density. Thus, we find that the current-current interaction in Galilean invariant 
theory is 

£>l = — X (j - jimp — f)]4imp ~ J4pimp) ■ (192) 

The corresponding Hamiltonian including both £7 and CJ I is 

d 3 x [-0 f ( — ^Jf)i> + -j^HH^] + A d 3 x pp lmp + X' / d 3 x (j ■ j iTOp - pj 4lmp - tePimp) • (193) 

The classical limit, where a pointdike impurity is moving along a given trajectory, is obtained by the following 
replacement : 

pim P {x, t) —> p c (x, t) = S(x - C(t)), (194) 

J imp (x, t) -> Jc (x, t) = C{t)6(x - C(t)), (195) 

kim P (x, t) -> j 4c (z, t) = ^CW 2 <5(a; - C(*))- (196) 

C. General solution for the current-current interaction 

The general solution within the Bogoliubov approximation is easily obtained for the Hamiltonian (193) in the same 
manner as before. After the Bogoliubov transformation, the Hamiltonian is 

H~E' Q + ^kb{{t)b k {t) + J2'(Fk(t)bi(t)+h.c). (197) 

Here Fk(t) is defined by 

F k (t)=nX. 



efc 



- u\ , A ' ^kk-jkit) X' e fe . ' 



(198) 



Fourier transforms of p c (x,t), j c (x,t), and j4 C (x,t) are denoted by f>k(t), Jk(t), and j4k{t), respectively, i.e., 

Pfe(i) = J d 3 x Pc (x,t)e- tk x , (199) 

J*(t) = / d 3 a;j c ( : r,t) e - ifea5 , (200) 

hk(t) = J d 3 xj4 C (x,t)e- ikx . (201) 

By solving the Heisenberg equation of motion, we obtain 

&L(*) = &£ + (*) + **(*). (202) 
&*(*) = &*(*) + **(*)■ (203) 

Here $&(£) is 

**(*) = ~e- iWkt [ dt' F k (t')e Wkt ' . (204) 
" J-00 

The same analyzes as were done before, will hold for the above general solution by replacing 4>k(t) with $k(t). 
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VIII. CONCLUSION AND DISCUSSION 

We have shown that the motion of classical impurities in the homogeneous BEC creates the time dependent coherent 
state. This results in macroscopic phenomena such as energy dissipation, induced density fluctuations, and so on. 
The essence of our results is that the homogeneous condensate is very robust against the motion of the impurities in 
the large N limit. However, in the current experiments where the typical number of atoms are on the order of 10 5 ~ 7 , 
the motion of the impurities could destroy the condensates by increasing the number of impurities. In order to study 
these effects, wc need numerical analysis with additional trapping potentials. An extension to the finite temperature 
is also desirable for this purpose. 

We have also shown that our simple model revealed various properties of the homogeneous condensate which 
have not been discussed before. These properties support the possibility of BECs being the space-time medium for 
impurities. Minkowski space-time structure was shown to be a certain part of the homogeneous BEC in the low 
energy scale, such as the wave equations (111), an effective Minkowski metric (104) and its inverse (105). In this 
context special relativity could be considered as an effective theory rather than a fundamental one, which emerges 
from the underlying microscopic degrees of freedom. Then it is quite natural to expect that the Lorentz symmetry 
is not the exact symmetry at all energy scales. We again emphasize the non-triviality of our result regarding the 
energy momentum stress complex. The deviation of this quantity from the homogeneous condensate is shown to be 
symmetric with an effective Minkowski metric (104). 

Several examples worked out in Sec. V illustrate the significance of our model. The effective attractive Yukawa type 
force between stationary impurities in our model is understood as the exchange of bogolons whose dispersion relation 
is linear. This effect itself seems difficult to observe directly with the current limitations in experimental techniques. 
Because this force decays exponentially for distance longer than the coherence length. The microscopic derivation of 
Landau's criterion shows that this criterion is not so strong one as was originally proposed. It is found that there 
occur finite amount of energy dissipation even slightly below the critical velocity given by eq. (119). However, this 
dissipation decays exponentially for these subluminal motion of impurities. We found that energy dissipation for the 
uniformly accelerated motion of impurity is too small to observe by direct measurements. In stead we suggest to look 
for a circular motion case in current experiments. Further investigations for this case need to be done with numerical 
analysis. 

There are several future directions in studying this problem. First, we need to treat the impurity at the quantum 
level rather than as a point-like classical object. As we have already seen, the dynamics of a single quantum point-like 
impurity is difficult to solve analytically. Another possible extension of our model is to consider a quantized field 
for the impurities as was outlined in the previous section. Then, we solve the field equations for massive bosons and 
the impurities self-consistently. Another issue which was not studied in our model is inhomogeneous condensates. 
Particularly, impurity effects in the presence of quantized vortices is of interest. We would like to continue to study 
these problems in near future. 

We finally remark that our result does not agree with the conventional approach to quantum field theory in curved 
space-time [32]. According to this work, a uniformly accelerated particle detector in the empty Minkowski space- 
time is shown to detect a thermal spectrum, which leads to an interpretation of the thermal effect. If we follow the 
idea proposed in refs. [7-15] and if there exists the Unruh effect, then the spectrum of bogolon for the uniformly 
accelerated impurity case should be the thermal one. As we have shown in our model, the bogolon spectrum for the 
uniformly accelerated impurity in the homogeneous condensate is not the thermal spectrum (eq. (138) or (140)). We 
note that similar, but phenomenological, considerations were presented in Chapter 6 of ref. [19] [56]. The author of 
this book claims that the excitations caused by the uniformly accelerated relativistic motion of a point-like classical 
impurity in a supcrfluid 3 He are emitted with the thermal spectrum. His argument is solely based on the simple 
Boltzmann factor evaluated in the WKB approximtion. The presence of this Bolzmann factor, not a Bosc-Einstein 
factor, for the transition rate leads that author to the conclusion that an analog of the Unruh effect is present in the 
B-phase of 3 He. However, we would disagree with the interpretation that the homogeneous condensate would serve as 
a thermal bath for accelerated impurities. The reason is clear as seen below [57]. The thermal nature of reservoir is 
a consequence of a statistical ensemble average. From this averaging procedure, the expectation value of occupation 
number for bosons results in the Plankian distribution, but not the other way around. In other words, an apparent 
mathematical expression does not mean that the system is in a thermal equilibrium at certain temperature. This 
point was also discussed by Lee within the context of the Unruh effect [57]. In our model the homogeneous condensate 
is seen as a coherent state by the bogolon, i.e. eq. (63). Hence, the occupation number (64) with respect to the 
homogeneous condensate has nothing to do with any statistical ensemble average. This observation leads us to the 
main conclusion that the homogeneous condensate does not serve as a thermal bath for accelerated impurities. If the 
homogeneous condensate is the physically correct model for Minkowski space-time, it then logically follows that the 
apparent thermal response of the linearly accelerated detector models may be the result of improper regularization. 
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APPENDIX A: RELATIVISTIC CIRCULAR MOTION 



We look for a relativistic circular motion in (2 + 1) space-time. We assume that the four acceleration w fJ -(s) = 
d 2 x^(s)/ds 2 ((1 = 0,1,2) is of the form: 

w°(a) = 0, Z(a) = w^s) + iw 2 {s) = %e^ s \ (Al) 

where s is the proper time, a is an acceleration constant, and ip(s) is a real function of the proper time. Thus 
w^Wfj, = —a 2 /c 4 and is a space-like. The four velocity u M (s) = dx fl (s)/ds satisfies u^u^ = 1. We solve this 
differential equation with the following initial condition, 



dx°(s) 

Z(a)\ a=0 = R , Z(s)\ s=0 = iR a n 1 ' 



ds 



(A2) 

s=0 



with given radius i?o and angular velocity Slo corresponding to the non-rclativistic limit. The relation between the 
acceleration a and these parameters Rq, flo is a = RqUq. The solution is 

x°(s) = ct(s) = 7s, (A3) 

Z{s) = x(s) + iy(s) = i?o7 2 cxp(«^) - R j 2 /3 2 , (A4) 

Rq7 



where (3 = Roflo/c and 7=1/ \J\ — (3 2 the relativistic factor. From this solution we find the trajectory of relativistic 
circular motion in terms of time t as 

x(t) = i? 07 2 [cos(^) - (3\ y(t) = R al 2 sin(^). (A5) 

In the limit c — > oo, this solution becomes 

x(t) = Rq cos(Q, t), y(t) = Rq sin(Qoi)- (A6) 
We remark the center of the circular motion is not fixed at the origin before the limit c — > oo. 
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